Resonant inelastic X-ray scattering in a Mott insulator 
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We calculate the resonant inelastic X-ray scattering (RIXS) response in a Mott insulator which 
is described by the Falicov-Kimball model. The model can be solved exactly within the single site 
dynamical mean-field theory (DMFT) approximation and the calculated RIXS response is accurate 
up to a local background correction. We find that on resonance the RIXS response is greatly 
enhanced over various other non-resonant background efltects and the response systematically evolves 
from a single peak structure, arising due to relaxation processes within the lower Hubbard band, to 
a two peak structure, arising due to relaxation processes within the upper Hubbard band as well 
as across the Mott gap into the lower Hubbard band, as we vary the incident photon frequency 
to access states from the bottom of the lower Hubbard band to the top of the upper Hubbard 
band. The charge transfer excitations are found to disperse monotonically outwards (as a function 
of transfered energy) as we go from the center of the Brillouin zone towards the zone corner. These 
correlation induced features have been observed by Hasan et. al. (Science 288, 1811 (2000)) and 
many other experimentalists in RIXS measurements over various transition metal oxide compounds 
and are found to be robust and survive even for large Auger lifetime broadening effects. As a 
comparison, we also calculate the dynamic structure factor for this model, which is proportional to 
the nonresonant part of the response, and does not show these specific signatures. 

PACS numbers: 71.10.Fd, 71.27.+a, 74.72.-h, 78.70.Ck 



I. INTRODUCTION 

Resonant inelastic X-ray scattering is essentially a deep 
core level spectroscopic method which is increasingly be- 
coming an essential technique in understanding the com- 
plex electronic dynamics of a wide class of novel materi- 
als like cuprates, manganites and various other transition 
metal oxide compounds. In the RIXS process, a highly 
energetic X-ray photon (with energy ~ 1 — 10 keV) ex- 
cites a deep core level electron into the unoccupied states 
of the conduction band. The excited electron then under- 
goes inelastic scattering processes with various intrinsic 
excitations present in the system and finally, a conduc- 
tion band electron fills up the core-hole and emits a pho- 
ton with relatively lower energy. So, this is a two photon 
inelastic process with no core-hole present in the final 
state. The transfered energy and momentum to the in- 
trinsic excitations of the system as well as the change 
in the polarization of the scattered photon can provide 
important information regarding these excitations. Also, 
RIXS being a resonant technique, the incident photon en- 
ergy can be chosen to coincide with, and hence resonate 
with, certain intrinsic X-ray atomic transitions which in 
effect can greatly enhance the inelastic scattering cross 
section and hence can be used as a probe for charge, mag- 
netic and orbital degrees of freedom on selective atomic 
sites. 

RIXS has several advantages over other spectro- 
scopic techniques like angle resolved photo-emission spec- 
troscopy (ARPES) and neutron scattering. First, in 
ARPES, the incident photon knocks out an electron from 
the system and hence can only probe the occupied states 
in a system but RIXS, being a high-energy process can 



excite a system into unoccupied intermediate states (like 
the upper Hubbard band in a Mott insulator) and hence 
can be used as a probe for understanding complex elec- 
tron dynamics in those strongly correlated intermedi- 
ate states. Inverse photo-emission techniques, in which 
an electron is injected into the system, can also access 
the unoccupied states of a system, but this method will 
charge the system and so far no momentum resolved in- 
verse ARPES with sufficient energy resolution has been 
developed. Secondly, the scattering phase space, i.e. the 
range of energies and momenta that can be transfered 
in the RIXS process, is much larger than other available 
photon scattering techniques involving visible or infrared 
light. As a result, RIXS can probe low-energy excita- 
tions over a wider range of the Brillouin zone and most 
importantly it can be used to probe all three directions 
of the Brillouin zone, and hence can be used even for 
materials which are intrinsically three dimensional in na- 
ture as compared to the ARPES which, because of the in 
plane momentum conservation, is widely used for mate- 
rials which are inherently two dimensional in nature and 
for three dimensional materials analysis of ARPES spec- 
tra are much more complicated (the perpendicular com- 
ponent of momentum is integrated over in the ARPES 
spectra). Also, photons can much more strongly couple 
with the electronic system through Coulomb interaction 
as compared to the neutron-electron interaction arising 
through the spin exchange interaction. As a result, RIXS 
can be used on very small volume samples, thin films, sur- 
faces and nano-objects in addition to bulk single crystal 
or powder samples. Also, for materials for which growing 
large volume single crystals is a difficult task, RIXS has 
significant advantage over neutron scattering. Besides 
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these, RIXS is polarization dependent and hence can be 
used to probe magnetic excitations and also it can be 
used as a probe for certain specific elements or orbitals 
in a system. The main disadvantage of RIXS over other 
techniques is, it requires substantially large incident pho- 
ton flux in order to have comparable or better energy and 
momentum resolution. But recent progress in RIXS in- 
strumentation has dramatically improved upon this sit- 
uation and RIXS is beginning to become an important 
probe for condensed matter physics. 

Over the past decade or so, RIXS has been carried out 
over large classes of transition metal oxide compounds 
like cuprates, manganites, irridates etc. (see Ref. fll 
for a detailed review). Most notably, RIXS measure- 
ments have been performed at the Cu K-edge^"— over 
a large class of cuprates like the undoped hd^^GnO 4^^^ 
Nd2Cu04,2''* SraCuOaClsr^ CazCuOzClz,^^ quasi one 
dimensional cuprates like SrCu02)^i^ Sr2Cu03}S, hole- 
doped cuprates La2-2;Sr2;Cu04,— YBa2Cu307_5,— and 
electron doped cuprates Ndi.85Ceo.i5Cu04fi^ as well 
as at the Cu La-edge—"— over various undoped 
cuprates CnOP Sr2Cu02Cl2^ La2Cu04,12J4. and 
doped systems Bi2Sr2CaCu208+5, Nd2-2:Ce2:Cu04,~ 
La2-xSr2,Cu04ji2r— . Besides these, RIXS measurements 
at the Mn K-edge in the orbitally ordered manganite 
LaMnOsfi^ at the Mn L2,3 absorption edge in MnO^^ 
at the Ni Ls-edge in NiOii and at the Ir Ls-edge in 
the bd Mott insulator Sr2lr0 4^^i^^ have also been per- 
formed. A common feature of these materials is that all 
of them are either Mott insulators or doped Mott insula- 
tors and have interesting magnetic ground states. RIXS 
measurements on these materials have probed energy and 
momentum resolved features of charge transfer excita- 
tions^i^ii^, dd excitationsi^ii^ii^ (arising due to transi- 
tions between crystal field split d orbitals), orbitons^^ in 
orbitally ordered systems and even magnetic excitations 
like magnons^ and bi-magnonsii. 

Theoretical approaches in understanding the RIXS re- 
sponse are mainly either based on exact diagonaliza- 
tion of model Hamiltonians over finite but small clus- 
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or based on a single particle approach^^i2i 
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ters 

which includes realistic band structure effects but treats 
the correlation effects perturbatively under the random 
phase approximation^^"— (RPA) or under a self consis- 
tent renormalization^S (SCR) approach and the effect of 
scattering from the core-hole in the first order Born ap- 
proximation or multiple scattering approximation. The 
exact diagonalization method treats the strong corre- 
lation effects exactly but because of the exponentially 
growing basis problem this method is limited to small 
size clusters and small number of orbitals and hence has 
limited momentum resolution. Also, the effects of the 
core-hole in this approach as well as in the SCR based 
calculations are taken either through an input core-hole 
lifetime, arising due to Auger and fluorescence effects, 
which broadens the intermediate states or under the 
ultrashort core-hole lifetime (UCL) approximation^ii^ 
which is found to be perturbatively exact for small as well 



as large core-hole potentials. But the effect of the core- 
hole lifetime arising solely due to intrinsic strong corre- 
lation effects in a Mott insulator on the RIXS response 
has not been addressed so far and in this work we use the 
Falicov-KimbalP'^ (FK) model to address this issue. The 
main motivation in choosing the FK model is that the FK 
model is one of the simplest models of strongly correlated 
electron systems which can be exactly solved^ii^^ under 
the single site dynamical mean field theory^^i^ (DMFT) 
approximation and most notably shows a Mott insulating 
ground state for large interaction strength between the 
itinerant and the static electrons. Also, the fully renor- 
malized two particle dynamic charge correlation function 
involving the itinerant species as well as the finite tem- 
perature core-hole propagator in this model can be calcu- 
lated exactly. In the following sections, we show that the 
calculated RIXS response in the limit of large core-hole 
energy is "exact" provided we : (1) neglect some mo- 
mentum independent background contributions and (2) 
calculate the charge vertex for exchange processes (see 
below) under the Hartree-Fock approximation, which is 
correct in leading order. We believe this model specific 
calculation can shed some light onto our understanding 
of the strong correlation effects in the RIXS response and 
the charge dynamics in a Mott insulator in general. The 
organization of the paper is as follows. In Sec. II, we pro- 
vide a brief mathematical formulation for the calculation 
of the RIXS cross section followed by Sec. Ill where we 
show a more detailed calculation for the RIXS response 
in the FK model. In Sec. IV, we show our results for a 
half filled Mott insulator followed by Sec. V where we 
discuss the core-hole lifetime broadening effects on the 
RIXS response and in Sec. VI we show some results for 
the case of particle-hole asymmetric Mott insulator. Fi- 
nally, in Sec. VII we conclude. 



II. MATHEMATICAL FORMULATION OF RIXS 

Our starting point is the familiar electron-photon in- 
teraction Hamiltoniai*^ 
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for a system of TV electrons. A(r, t) is the vector po- 
tential for the external electromagnetic field and can be 
expanded in a plane wave basis as 



A(r, i) = 51 VA^ (eak.e'^''-'-'^*) + e*al^^t 



,(2) 



where An — tttj^ — . V is the volume of the system, e is 
the polarization of the light and we have fixed the gauge 
by choosing V • A(r, i) = in Eq. (P). 
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In the RIXS process, an incident X-ray photon with 
momentum , energy uji and polarization is scattered 
to a final state described by momentum k^, energy w/ 
and polarization Sf. Fermi's golden rule to second order 
in Hint gives the transition rate for this process: 
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where |I) = \i) (g) \ki,uji), |F) = |/) » |k/,w/) and |7i) 
are the matrix product states for the initial, final and 
the intermediate states of the systems respectively (both 
electronic and photon states are present in the initial and 
final states while the intermediate states are just the elec- 
tronic states) and Ei, Ep and e„ are the corresponding 
energies, respectively. It is interesting to mention that 
the intermediate state |n) has a core-hole while the initial 
and final states, |I) and |F) have no core-hole. The first- 
order amplitude is in general dominant over the second- 
order contribution except near resonance when the inci- 
dent photon energy is nearly equal to a specific atomic 
transition in a material; i. e., coi « e„ — e^. At reso- 
nance, the second-order term becomes overwhelmingly 
large compared to the first-order term and hence the 
second-order term causes resonant scattering while the 
first order term gives rise to nonresonant scattering. 

The diamagnetic term proportional to as well as the 
spin-orbit coupling term proportional to cr • (dA/dt) x 
A in Eq. ([T]) contribute to the first-order amplitude. 
The latter is smaller than the former by a factor of 
u!i^ff/mc^ ^ 1 and also, at resonance, the contribution 
from the diamagnetic term is negligibly small compared 
to the resonant term and hence their contribution will be 
neglected. So, then the resonant part of the second-order 
amplitude at zero temperature is given by^ ''^^'^^ 
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where = — cjj is the transfered energy and 
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(5) 



is the relevant transition operator for the RIXS cross sec- 
tion. The first term in Eq. ([5]) causes nonmagnetic scat- 
tering while the second term, arising from the spin-orbit 
coupling term in Hint , causes magnetic scattering which, 
for typical incident photon energy ('^ 1 — lOkeV) and the 
localized core levels involved in a RIXS process, is about 
100 times smaller than the non-magnetic term^ and hence 
will also be neglected. Finally, under such circumstances, 
we assume the dipole limit for the RIXS process and the 
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FIG. 1: The Kadanoff-Baym-Keldysh contour. The contour 
starts at time 0, moves forward along the real axis to time t 
then moves backward along the real axis to time and finally 
downwards along the imaginary axis to time —ijS. 



transition operator is then given by 



V = e 



1 ^ 

■ D with D = V p. 



(6) 



III. RIXS RESPONSE IN THE 
FALICOV-KIMBALL MODEL 

The single-site Hamiltonian of the Falicov-Kimball 
model (in the hole representation) including the inter- 
action with a core-hole is given by 

Hioc = UudUf + QdUdUh + QfUfUh - ^J-nd 

+ {Ef~fi)nf + {Eh~fi)nh (7) 

where Ud — Sd, nf = f^f and Uh = h^h are the oc- 
cupation number operators for the d-hole, /-hole and 
core-hole state, respectively. U is the onsite Coulomb 
interaction between the itinerant d and the static / 
holes, Qd > and Qf > are the Coulomb interac- 
tions between the core-hole and the d-hole and /-hole, 
respectively, Ef is the site energy of the / state and 
Eh ^ 10^ — 10* eV is the energy of the core-hole state. 

The full Hamiltonian on the lattice includes a repeat of 
this local Hamiltonian for each lattice site and a hopping 
of the itinerant d-holes between neighboring sites. The 
density matrix for the single impurity problem in DMFT 
is then given by 



rcexp<^ / dt' dt"d''{t')Mt',t")d{t") \ ,(8) 



P 



where the time-ordering and integration are performed 
over the Kadanoff-Baym-Keldysh contour as shown in 
Fig. [T] and /3 = l/fc^T is the inverse temperature. Here, 
the dynamical mean- field Xc{t',t") and the chemical po- 
tential fi are taken from the equilibrium solution of the 
conduction electron problem without the core-hole, aris- 
ing under the single site DMFT approximation. This 
in effect implies that we are treating the creation of the 
core-hole under the sudden approximation instead of a 





FIG. 2: (Color online) RIXS interaction vertices. Wavy lines 
(blue) represent incident or scattered photons, the dashed 
lines (green) represent the core-hole propagator, and the solid 
lines (red) represent the propagator for the itinerant d-holes. 
The labels indicate the energies of the different particles. Note 
that, the momentum and energy are conserved at each vertex. 



full non-equilibrium treatment of the core-hole propaga- 
tor. 

The impurity problem arising under the DMFT ap- 
proximation can be solved exactly in this case and the 
local d-hole propagator, G'^''{uj), is given by 



Wq 



Wi 



,(9) 



where uj^ — u) US {6 > 0) , Wq and Wi are the probabil- 
ities for finding a given site unoccupied and occupied by 
an /-hole, respectively. The momentum dependent fully 
renormalized d-hole propagator is given by 



Gd(q,w) = 



1 



(10) 



where the local self energy is related to the local propa- 
gator through Dyson's equation 

S^-(..+) =u;++f,- A(c.+) - [Criu^^)] ' ■ (11) 

Similarly, the core-hole Green's functions, G^{t) = 
-i{h{t)h^O)) and G^{t) = i{h{t)h'f {0)) can also be cal- 
culated^ by using either numerical integration over the 
Kadanoff-Byam-Keldysh contour or by the Wiener-Hopf 
sum equation'^^''"^ approach [the angular brackets ( ) de- 
note a trace over all states weighted by density matrix 
in Eq. [S] and the operators are in the interaction repre- 
sentation with respect to Hiod- Also, it is important to 
mention that, for the calculation of the itinerant as well 
as the core-hole propagators we use the d dimensional 
hypercubic lattice density of states (DOS) in the limit of 
d — oo (DMFT approximation) . 

The interaction of the X-ray photon with the electronic 
subsystem of matter can be represented by the diagrams 
shown in Fig. [2l We have explicitly shown the direct 
dependence of the core-hole propagator on the core-hole 
energy Eh, which is supposed to be much larger than the 
band energies and is of the order of the incident (wi) and 
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■UJf 
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Gdd(q,ix;) 
- *- - Gu,{uj) 

FIG. 3: (Color online) Bare loop contribution to the direct 
RIXS process. The bottom diagram gives a momentum- 
independent contribution (see text) to the RIXS process and 
hence will be neglected in this calculation. 



the scattered X-ray photon energies (w/), respectively. 
One can see that in the case of large photon and core- 
hole energies only the first two vertices (top two diagrams 
in Fig. contribute significantly whereas the contribu- 
tion from the remaining two (bottom two diagrams in 
Fig. [2]) are negligibly small because the hole-propagators 
are evaluated too far off the energy shell. 

In that limit, the bare-loop contribution to the am- 
plitude for the RIXS process is represented by the two 
diagrams shown in Fig. [3] The contribution of the 
top diagram to the four-particle correlation function, 

f ii^i^ii ii^f, ii^'f, evaluated on the imaginary axis, 
is equal to 

- \ '^xt'^{i^m - iVi + iVf,iuJra\^) 

m 

xxQ\iu:yn + ivf,iu:ni + w'f) (12) 

whereas the contribution of the bottom diagram to 
X^ljj,i{wi,ivf,iv'j,iv'^) is equal to 

1 V ^ hh / ■ ■ ■ \ 

771 

xxo''(iw,„ - ii//,ia;„i - w}|0), (13) 
where we have introduced the bare charge susceptibilities 



IV, 



xGd(k,iu;„), (14) 
X^o'^iiu^t + ^^t) = -Ghiiojt + ^y)Gh{^u^t) (15) 
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and iw^ = iuJm — Eh- Since the core- hole propagator is 
local, the bottom diagram in Fig. [3] does not depend on 
the photon wavevector and hence can only contribute to 
momentum-independent background effects as evident in 
Eq. (jl3l) , whereas the first diagram depends on the trans- 
fered momentum q = ki — k/ . As in the present study we 
are interested in the energy and wavevector dependence 
of the RIXS response, we neglect all such momentum- 
independent contributions. The technical reason behind 
neglecting these terms is that the momentum indepen- 
dent contributions are the local ones and hence they in- 
clude all types of many body scattering processes and can 
only be derived from the solution of the single-impurity 
problem for the four-particle correlation function x^^^ 
(which involves multiparticle vertices and many more 
complications) and at this moment we do not have well 
developed approach for this. 

It is important to mention that RIXS processes^ii^ 
can happen either through a direct process, as in the case 
of the L2,3-edge 2p — > 3d RIXS, in which the core-electron 
is excited to an unoccupied state of the correlated valence 
band (d-band) or through an indirect process, as in the 
case of K-edge Is 4p RIXS, in which the excited core- 
electron goes into an uncorrelated Ap band several eV 
above the Fermi level. In the present case, the excited 
core-electron goes to the correlated d-band and dipole 
selective transition 2p 3d or 3p — > 4d is consistent 
with the involvement of a 2p or 3p core-hole and hence 
our study is related to the direct RIXS processes like L- 
edge or M-edge RIXS. In particular, we will be studying 
the L-edge RIXS in the following sections. 

The calculation of the RIXS response involves analytic 
continuation to real frequencies, which is a well defined 
but tedious procedure'*^ and gives the following contri- 
bution to the RIXS cross section, 



1 

'2^ 



+ 00 



du: [f{Lj) - f(uj + Q)] xT{<^- + c^,;, + w.) 



xRelxi^iu-,^- +^\ci)-xi''{^+,^- , (16) 
where /(w) = l/[exp(/3w) + 1] is the Fermi function and 

x'$'\lO- +LO^,LO+ = - \Gh{LO+ + L0^ - E^f . (17) 

Similary, the contribution from the second diagram is 
given by 

— / [fiiu) - fiuj + n)] xi\oJ- +uJf,uJ++ c^.lO) 

xRe [xTiio- ^n,Lu-)^ Xo"(^+ - ^,^')] (18) 

and as already has been mentioned it does not depend 
on the photon momentum and hence will be neglected. 

Next, we introduce the renormalization of the bare 
charge susceptibilities through inclusion of charge ver- 
tices. In the simplest case, this can be done by insert- 
ing the two-particle charge vertex either in between the 
two d-hole propagators (dd-channel) or between the two 




FIG. 4: (Color online) Diagrams for the RIXS response. Panel 
(a) total direct resonant scattering, (b) full exchange resonant 
scattering and (c) partial exchange resonant scattering pro- 
cesses. Explicit Eh dependence of the core-hole propagator is 
not shown in the figure. 



core- hole propagators (/i/i-channel). First, we consider 
the effect of charge screening in the dd-channel which is 
shown in the second diagram in Fig. 2] (a). The sum 
of the two diagrams in Fig. SJa) corresponds to the re- 
placement in Eq. ([16]) of the bare charge susceptibility 
Xo'^{ioJm + ii^,i(^m\(l) by the fully renormalized charge 
susceptibility x'^'^i'i'^m + i^^, iwm|q) which in the case of 
Falicov-Kimball model is given by 

[X'^'^iil^m + it^mlq)] ^ = [xo''(it^m + «^^m|q)] ^ 

+T{iLOm.+iv,iUJm)-, (19) 

where the irreducible charge vertes^*^ is 



r(jlLJ,„ + iv, iuJrn) 



1 S(icjm) — E(?w„i + iu) 

T G{iuJm) - G{iuJrn + w) 



(20) 



The sum of the two diagrams in Fig. IH^a) corresponds 
to the direct scattering contribution to the RIXS process 
and is given by 

1 /■+°° 9 
--^ / dLo[f{Lo)- f{LO + n)]\Gh{u:+ +u:,-Eh)\ 

xRe[x''\u:-,u:- +n\ci)-x''\u:+,,o- +^\ci)\ . (21) 

The word direct here refers to a direct Coulomb scat- 
tering process since we already discussed that we will 
calculate only direct RIXS process in this work. The 
contribution of the diagram in Fig. IDJb), which we la- 
bel as the full exchange scattering processes to the four- 
particle correlation function, x^ff / ^'^/i ^^^i); on 
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the imaginary axis is given by 



xx'*'*(^l^,;-^^^/,q), (22) 



where 

dd 



27r 



+c>o 



^X^'i^ - ^+|q) - X''(^ - ^-|q)] /(^), (23) 
is the dynamical charge susceptibihty and A satisfies 

K{ivi,iuf,ivi-ivf) = - '^'^Gh{i^^m + i^i - E^) 

m 

xGdiiuj„i)Gh{iuj,n + ivf - £^/i).(24) 

It is important to mention that in the derivation of the 
full exchange scattering process contribution to RIXS 
in Eq. (1^^ . we have approximated the core- hole-d- hole 
charge vertex by the core-hole - d-hole interaction Qd 
under the Hartree-Fock approximation, which is exact to 
leading order. After analytic continuation to real fre- 
quencies, the corresponding contribution to the RIXS 
cross section is equal to 

— Im [k{uji + i5,ujf + i6,VL + iS) 



xA{uji - iS,ujf ~ iS,n + i6)x'^'^{^ + iS, q)] 
where A(a;i ± oj^ ± i(5, 57 -I- i5) is given by 



(25) 



+ 00 



du f{Lu) [Gniio^ + Lu1)Gh{L^^ + ij'})lmGd{Lo+) 

+Gdiuj^ - u;^)Gh{uJ^ -n- Eh)lmGhioj+ - En) 
+Gd{^^ - ujf)Gh{^^ EH)lTnGh{<^+ - Eh)] (26) 

and Lu^ ^ = LUij — Eh are the incident and scattered pho- 
ton energies measured with respect to the core-hole en- 
ergy. Eh- 

For large core-hole energy. Eh, we need to keep only 
the first term which has a small difference of energies 
fi + ujij — Eh and can safely neglect the other two terms 
containing a large differences in energies {fi — Eh and 4- 
il — Eh, respectively). We then obtain the contribution 
from the full exchange processes as 

\A{uj, + i6,ujf + iS)\ liRx'^'^in + i6\q) (27) 



where. 



+ 00 



(28) 



A{uji + iS,uj f + iS) = — doj f {uj)G h{i-^J^ 

xG,i(a;+ -f a;^)ImGd(cj+), 
lmx''\n + tS,q) = — J di,[f{Lu)-fiuj + n)] 
xRe [x''^{ij-,ij- + n\cD - x'^'^(w+,w" + 57|q)] . (29) 



Apart from the direct and full exchange resonant scat- 
tering processes, we also have processes, termed as the 
partial exchange processes, as shown in Fig. |4ljc) [this 
process arises due to the cross terms in the square in 
Eq. ([3])]. The sum of their contributions to the four- 
particle correlation function, ^ ^(w^, w/, w^, if^'), on 
the imaginary axis, is equal to 

- Qd-n'^x''''''i'i-^m7'i-^m + ivi - ws\oi) 

X [GhiiuJm + iv'i - Eh)A(ivi,ivf, ivi - ivf) 
+Gh{ii^ni + w% - Eh)A(iv[, iv'p iv[ - iv'j)\ . (30) 

After analytic continuation to real frequencies, we obtain 
the following partial exchange contribution to the RIXS 
response 

~^Rei^A{uj,+iS,LUf + iS) J duj[f{uj)~ f{uj + n)] 
xGhiuJ- + - Eh)Re [x'^'^ {cu- , lj- + n\q) 



-X'^uo' 



,L0- +n\q)\]. (31) 

It is interesting to mention that for the i-edge RIXS 
process, which we consider in the present study, pure 
direct resonant scattering processes are overwhelmingly 
dominant over the resonant exchange (both full and par- 
tial) processes. It is also interesting to mention that in 
Fig. 4(a) inclusion of the reducible charge vertex, which 
in general is nonlocal, into the core-hole-core-hole chan- 
nel will give a vanishing contribution to the RIXS cross- 
section. This is because of the local nature (momentum 
independence) of the core-hole propagator and the van- 
ishing of the dd-bubble (charge susceptibility) in the dia- 
grams like the one presented in Fig.Ufb), but aligned (ro- 
tated) in the vertical direction (with the external photon 
lines kept unchanged), in the uniform limit (q = 0). Be- 
sides the diagrams considered above, we could also have 
considered other contributions to the charge vertex be- 
tween the two core-hole propagators which produce dia- 
grams like the parquet diagram in Fig. [5] From a simple 
power counting argument we can show that the contri- 
bution of such diagrams goes at least as an inverse power 
of the dimension of the lattice, d. So, in the limit of 
d ^ oo, they all have vanishing contributions except in 
the case when they are all local (and we neglect all such 
momentum- independent contributions). 

We could also have included renormalization through 
the core-hole loops 



+ 00 



dujf{Lu)lmGh{i^+ ~ Eh) 



X [Gh{uj + iv- Eh) + Gh{uj -IV- Eh)] (32) 

but they also contain large differences of energies fi — Eh 
and fi + fl — Eh and hence can be neglected for large 
core- hole energy. Eh- 
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FIG. 5: (Color online) Parquet diagrams with vanishing con- 
tribution to the RIXS response in the limit of infinite dimen- 
sion. 

The nonresonant part of the RIXS response is found 
to be related to the density-density correlation function. 
To be precise, the nonresonant part is proportional to the 
dynamical structure factor^, S{q, fi) which is given by 

Siq, n) = ^-[l + nein, T)] Irnx'^^in + i6\q), (33) 

TT 

where nB{oj,T) = l/[exp{Puj) — 1] is the Bose distribu- 
tion function and x'^'^i^ + *<5|q); is the dynamic charge 
susceptibility of the system. 



IV. RIXS RESPONSE FOR THE HALF FILLED 
MOTT INSULATOR 

As has been already stated, the Falicov-Kimball model 
at half filling {rif = 0.5) shows a Mott insulating ground 
state for ?7 > C/^ = v^. We choose U = 2.0, Qj^ = Qf ^ 
2.5 and T — 0.1 in units of effective hopping amplitude 
t* and from here onwards we choose t* = 1. This choice 
of U gives a small "gap" (Agap — 0.25 in units of t*) 
Mott insulator. It is interesting to mention that for the 
d — >■ c» hypercubic lattice DOS there is no true gap as 
there is an exponentially small DOS inside the gap. 

In Fig. [6l we show the spectral function for the itiner- 
ant species, Ad{u}) as weU as for the core-hole, Ah{'jo — 
Eh). Ad{oj) clearly shows a "gap" (Mott gap) at the 
Fermi level (w = 0) while Ah{uj — Eh) also shows a "gap" 
at some other frequency and the origin of this gap is re- 
lated to the same strong correlation effects that gives rise 
to the Mott gap in the itinerant species spectral function. 
Surprisingly, the gap structure is quite different for the 
core-hole, which is arising due to the asymmetry in the 
Green's function for large Eh- Ah{ui) is dominated by a 
broad feature, arising from the projection of G^{t) onto 
the Uh = Q,nf = Q configuration, along with a very sharp 
peak arising from the projection onto the = 0, ti/ = 1 
configuration in the final state. 

From the knowledge of the itinerant electron prop- 
agator, G'd(q, w), core-hole propagator, G/i(w) and 
the fully renormalized two-particle charge susceptibility, 
X'^'^iiuJm, i^m + 'ii^lq), we can calculate the RIXS response 
either as a function of transfered energy (fi) for a given 
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FIG. 6: (Color online) Spectral function for the d-hole, Ad{io), 
and the core-hole, ^^(a;). For ^^(aj), the frequency is mea- 
sured with respect to the core- hole energy. Eh, as shown in 
the top horizontal axis. 

fixed incident photon energy, Ui (measured with respect 
to the core-hole energy Eh), or as a function of uji for 
a given fixed transfered energy, 51, for various transfered 
momenta q of the photon. It is interesting to mention 
that in the limit d ^ oo the momentum on the hypercu- 
bic lattice only enters through the dimensionless param- 
eteri^, 

1 

X(q)= hm -^cos(g.)- (34) 

i—1 

So, -1 < X < 1 and X = 1 and X = -1 corre- 
sponds to the center, (0, • • • ,0) (F point), and the cor- 
ner, (tt, • • • , tt) (M point), of the Brillouin zone of a d- 
dimensional hypercubic lattice, respectively. It is conve- 
nient to think of this parametrization as corresponding 
to RIXS scattering in the diagonal (1 • • • 1) direction. 

First, in Fig. [71 we plot the resonant part of the RIXS 
response as a function of the transfered energy, 57, for 
various transfered momenta, X, for three different inci- 
dent photon energies, uji = —0.5, 0.5, 1.5. For uJi = —0.5 
the incident photon energy is lying in the lower Hubbard 
band and the inelastic relaxation processes can happen 
only within the lower Hubbard band which is evident 
in the single peak structure in the RIXS response in 
Fig. El^a). At the M point the peak is large and well 
defined but as we go towards the middle of the Brillouin 
zone {X = 0.0) the peak gets broadened and the position 
of the peak does not disperse significantly. Finally, as we 
approach the zone center the position of the peak dis- 
perses significantly and moves towards lower energy, also 
at the same time the peak gets more and more well de- 
fined though the integrated spectral intensity under the 
peak gradually diminishes and eventually goes to zero at 
the center of the Brillouin zone {X = 1) which is related 
to the vanishing of the uniform charge susceptibility. 

As we increase the incident photon energy above the 
Mott gap, we start to excite the system into the upper 
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FIG. 7: (Color online) RIXS response as a function of trans- 
feree! energy (0.) for various fixed transferee! momenta (X) 
witli three different incident photon energies, u;. In pane! (a), 
the incident photon energy is at the lower Hubbard band, in 
(b) at the bottom of the upper Hubbard band and in (c) it is 
at the middle of the upper Hubbard band. The inset of each 
plot schematically depicts the possible RIXS relaxation pro- 
cesses. The chosen parameters are U = 2.0, Qd ~ Qf = 2.5 
and T = 0.1. Note, the signal vanishes exactly at X = 1. 



Hubbard band and in Fig. [Tfb) , we show a characteristic 
response when the incident photon energy Ui = 0.5 is at 
the bottom of the upper Hubbard band. Near the zone 
corner the response still shows a single peak structure al- 
beit shifted by the insulating gap, but as we go towards 
the zone center, the response at low energy develops a 
very narrow secondary peak separated from the broad 
main peak by the Mott gap. The low-energy peak arises 
due to relaxation processes within the upper Hubbard 



band and is non-dispersive in nature, whereas the high- 
energy peak arises due to relaxation processes across the 
Mott gap into the lower Hubbard band and is found to 
be dispersive over the Brillouin zone. With further in- 
crease in the incident photon energy to = 1.5 the low 
energy peak and the spectral weight under it, as shown 
in Fig.[71[c), grows significantly and is visible for all mo- 
menta along the (1 • • • 1) direction. However the intensity 
of this low-nergy peak shows non-monotonic behaviour - 
it first increases up to X = and then starts to de- 
crease and eventually vanishes at the zone center, X = 1. 
On the other hand the intensity of the high-energy peak 
monotonically decreases as well as disperses to lower en- 
ergies as we go from the zone corner towards the zone 
center and eventually vanishes at the center of the Bril- 
louin zone. 

In Figs.[8{a), (b) and (c), we show the systematic evo- 
lution of the RIXS response for three momenta X = —1, 
X = and X = 0.9, respectively. As we vary the inci- 
dent photon energy from the bottom of the lower Hub- 
bard band to the top of the upper Hubbard band the 
RIXS response evolves from a single-peak structure to a 
two-peak structure and when uji is beyond the edges of 
the Hubbard bands the response vanishes quickly due to 
an exponential reduction of the density of states which 
in effect drastically reduces the phase space for inelastic 
scattering. Also, the overall response decreases as we go 
towards zone center. 

In Fig. [HI we show a more detailed and systematic evo- 
lution of the RIXS response with varying incident energy 
in the upper Hubbard band. As we can clearly see, the 
low-energy peak which arises due to relaxation processes 
within the upper Hubbard band does not disperse over 
the Brillouin zone (except its intensity varies) whereas 
the high-energy peak (which arises due to relaxation pro- 
cesses across the Mott gap) shows significant dispersion 
over the Brillouin zone. This feature has been observed 
in RIXS measurements on a Mott insulator Ca2Cu02Cl2 
by Hassan et. al^ and has been attributed to strong 
correlation effects. 

It is interesting to mention that at zero temperature 
in the half-filled Mott-insulating ground state, the chem- 
ical potential as well as the Fermi level lies within the 
Mott gap which results in a completely filled lower Hub- 
bard band (LHB) and a completely empty upper Hub- 
bard band (UHB). As a result, when the incident photon 
energy is within the LHB the RIXS response will vanish 
due to the unavailability of any unoccupied state in the 
LHB to which the core electron can be excited and when 
the incident photon energy is within the upper Hubbard 
band the RIXS response will have a single peak structure 
arising due to particle-hole excitations across the gap in 
the presence of the core-hole potential. At finite tempera- 
ture some of the states near the top of the LHB thermally 
excite across the gap and occupy the bottom of the UHB. 
So, at finite temperature, if the incident photon energy 
is within the LHB, the core electron can be excited to 
thermally excited unoccupied states of the LHB, giving 
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FIG. 8: (Color online) RIXS response as a function of trans- 
feree! energy (Q) for various fixed incident photon energies, uui, 
varying from the bottom of the lower Hubbard band to the 
top of the upper Hubbard band for three different momenta, 
X. Panel (a) zone corner (M point) X — —1, (b) somewhere 
in the middle of the zone X = and (c) near the zone center 
X = 0.9. Inset of each plot shows the position of each X 
inside the first Brillouin zone. All other parameters are the 
same as in Fig. [T] 



rise to a peak corresponding to the relaxation processes 
within the LHB and if the incident photon energy is in the 
UHB then the core electron that is excited to the empty 
states of the UHB can undergo relaxation processes ei- 
ther through scattering from the thermally excited UHB 
electrons occupying the bottom of the UHB giving rise 
to a low-energy nondispersive peak or through particle- 
hole excitations across the gap giving rise to a dispersive 
high-energy peak. Also, the whole structure gets shifted 
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FIG. 9: (Color online) Systematic evolution of the RIXS re- 
sponse as a function of transfered energy (Q,) for six incident 
photon energies from the bottom of the upper Hubbard band, 
panel (a) uji — 0.5, to the top of the upper Hubbard band, 
panel (f) uji = 3.0. The dispersive nature of the high energy 
peak (Mott gap excitation) is noticeable. All other parame- 
ters are same as in Fig. [T] 



to higher energy with increasing incident photon energy 
due to the fact that with increasing incident photon en- 
ergy the transfered energy must also increase in order 
to have resonant scattering from the thermally excited 
states which predominantly occupy the bottom of the 
UHB and the top of the LHB. In the case of large gap 
Mott insulator, as will be shown in a following section 
(Sec. VI), the intraband relaxation processes from the 
thermally excited states are negligibly small compared 
to the excitations across the Mott gap. 

Finally, in Figs. [TOlJa), (b) and (c), we plot the RIXS 
response as a function of the incident photon energy, w^, 
for various fixed transfered photon energies, D, for three 
transfered momenta X — —1, X — and X — 0.9, re- 
spectively. For small ft, the RIXS response shows a two- 
peak structure in uJi which corresponds to the relaxation 
processes within the individual bands (upper and lower 
Hubbard bands). As we increase fl an additional peak 
develops between the two peaks. This peak corresponds 
to the interband relaxation processes across the Mott gap 
and grows very rapidly with increasing fl while the other 
two peaks decrease in intensity until we are finally left 
with a lone peak. Hence, we can infer that in a Mott 
insulator, interband relaxation processes across the Mott 
gap are dominant over intraband relaxation processes. 
Also, as we go from the zone corner to the zone center 
the intensity of the peaks decreases as in Fig. [71 

We also have calculated the dynamical structure fac- 
tor, S'(q, ri), which is proportional to the nonresonant 
part of the RIXS response^. In Fig. [ill we plot ^(q, n) 
for the small-gap insulator at T — 0.1. Near the zone 
corner, S'(q, 51) has a broad midinfrared peak but as we 



10 





FIG. 10: (Color online) RIXS response as a function of inci- 
dent energy, uJi, for various fixed transfered energies, Q,, for 
three characteristic transfered momentum, X. Panel (a) zone 
corner (M point) X = —1, (b) in the middle of the zone X = 
and (c) near the zone center X = 0.9. The inset of each plot 
shows the position of each X inside the first Brillouin zone. 
All other parameters are the same as in Fig. [T] 



go towards the zone center, a secondary peak develops 
near = and the integrated spectral weight under the 
midinfrared peak decreases and finally exactly at the zone 
center {X = 1) the midinfrared peak completely vanishes 
while the peak around = turns into a delta function 
which again arises due to the vanishing of the uniform 
charge susceptibility at finite frequency. This behavior 
is quite different from the resonant response where the 
two-peak structure is most prominent near the zone cor- 
ner and also the midinfrared peak in S'(q, O) is much 
smaller and much less dispersive than the high-energy 
peak observed in the resonant response and most impor- 
tantly the position of the peak cannot be identified with 
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FIG. 11: (Color online) Dynamical structure factor, ^(q, f2) 
for U=2.0 and T = 0.1. Note, it vanishes for X = I. 



any particular X-ray transition process. 



V. CORE-HOLE BROADENING EFFECTS 

In the preceding section, we have not included any ad- 
ditional core-hole lifetime broadening effects which can 
arise due to various nonradiative Auger and fluorescence 
effects and are important in the transition metal RIXS 
processes (we only included the intrinsic many-body ef- 
fects in determining the core- hole lifetime). In our cal- 
culation, we can easily include such effects by simply 
making the core-hole energy Eh complex i. e. by 
making the transformation Eh — ^ Eh — iV into the re- 
tarded Green's function, Gl{t) = Q{t)[Gl{t) - G^{t% 
and Eh — Eh + iT into the advanced Green's func- 
tion Gl{t) = e(-t) [G^{t) - G^{t)] , respectively. Here, 
r = - with T being the core-hole life time and for most 
materials F 100 - 400 meV. 
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FIG. 12: (Color online) Core-hole spectral function, y4h(a;), 
evolution with core-hole hole broadening parameter, P. All 
other parameters are same as in Fig. [T] 

First, in Fig. [121 we show the systematic evolution of 
the core-hole spectral function, ^^(w), with additional 
core- hole broadening effects parametrized by F. With in- 
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FIG. 13: (Color online) Evolution of RIXS response as a func- 
tion of the transferee! energy, il, and various fixed transfered 
momenta, X, for incident photon energy uji = 2.5 with var- 
ious core-hole lifetime broadening parameters, F. Panel (a), 
response for F = and the inset of panel (a) shows quasielastic 
peak in the blown up far-infrared region. Panels (b), (c) and 
(d) show responses for broadening parameters much smaller 
(F = 0.1), comparable (F = 0.2) and larger (F = 0.3) than the 
Mott insulating gap A^ap ~ 0.25. All other parameters are 
the same as in Fig. [T] For all the four sets, calculations have 
been performed from X = —1.0 (top most curve) to X = 0.9 
(bottom most curve) with a momentum interval 0.1. 



creasing F, the height of the sharp peak in Ah (w) reduces 
while its width significantly broadens but the asymmet- 
rical nature of the structure is largely retained. Also, the 
Mott gap in the F = 0.0 case (near uj — Eh ^ —1.0) is 
replaced by a dip in the spectral function and the tail of 
Ah{i^) increases with increasing F. 

In Figs. [12] (a), (b), (c) and (d), we show a detailed 
evolution of the RIXS response with various broadening 
parameters, F (measured in units oft*). We choose four 
characteristic parameters F — 0.0, F — 0.1, F = 0.2 and 
F = 0.3 corresponding to no broadening, much smaller, 
comparable and larger broadening compared to the in- 
trinsic Mott-insulating gap in the system, respectively. 
All other parameters chosen are the same as in the pre- 
vious case. In Fig. [T3] (a), we plot the RIXS response 
without any additional core-hole broadening effect and 
the response shows a clear two-peak structure with the 
peaks well separated by the Mott gap. The low-energy 
peak is nondispersive but the high-energy peak is highly 
dispersive in momenta and also changes its shape signif- 
icantly as we go from the zone corner to the zone center. 
Also, interestingly, close to the zone center a quasielastic 
peak develops in the far-infrared region. The intensity as 
well as sharpness of this peak increases as we go towards 
the zone center. 

In the presence of additional core-hole lifetime broad- 
ening effects, we no longer see a clear gap structure in 



Figs. [13] (b), (c) and (d) but the two-peak structure is still 
clearly evident for all momenta. Also the high-energy 
peak still remains dispersive throughout the Brillouin 
zone, while the low-energy midinfrared peak is more or 
less nondispersive in nature as in the case with no addi- 
tional core-hole broadening. The height of the quasielas- 
tic peak close to the zone center decreases with increas- 
ing F, which is related to the reduction of the overall 
RIXS response with decreasing core-hole lifetime, while 
its width and dispersive features remains similar to the 
F = 0.0 case. However, in a real experiment, finite resolu- 
tion and the resolution broadened tail of the huge elastic 
peak will mask such quasielastic features which will then 
not be observable. So, the presence of additional core- 
hole broadening effects can significantly modify the over- 
all response, but the most important qualitative features, 
like the two-peak structure and the dispersive features of 
these peaks remain similar. 



VI. RESPONSE AWAY FROM HALF-FILLING 

Finally, we consider a Mott insulator at arbitrary fill- 
ing to examine the breaking of particle-hole symmetry in 
the RIXS response. We choose U = 4.0, Qd = Qf = 5.0, 
Uf — 0.25 and T — 0.1. This choice of parameters gives a 
large-gap Mott insulator (Ag^p ~ 1.8 in units of t*). In 
Fig. [Wal. we show the d-hole spectral function, Ad{uj), 
as well as the core-hole spectral function , Ah{io), with 
(F = 1.0) and without (F = 0.0) core-hole broadening 
effects. Ad{uj) is dominated by two asymmetrical peaks 
separated by a large gap at the Fermi level. Ah{io) [with- 
out any additional core-hole broadening effect (F = 0.0)] 
shows two very closely spaced sharp peaks on top of a 
broad feature but with the inclusion of large core-hole 
broadening (F = 1.0) the whole structure gets drastically 
modified and Ah{id) resembles a broad nearly symmetri- 
cal single peak. 

In Fig. [HTb) , we show the RIXS response as a func- 
tion of the transfered energy, il, for various fixed trans- 
fered momenta, X, for F = 0.0. The response in this 
case is overwhelmingly dominated by a huge peak aris- 
ing due to relaxation processes across the Mott gap into 
the LHB, while the intraband (within the UHB) relax- 
ation processes as shown in the inset of Fig. [Hf b) are 
much weaker than the interband processes. As already 
has been mentioned in a previous section (Sec. IV), in 
a large-gap Mott insulator the density of thermally ex- 
cited states which occupy the bottom of the UHB is ex- 
tremely small and hence cannot provide any significant 
relaxation to the core electrons excited to the UHB. The 
high-energy peak, just as in the case of particle-hole sym- 
metric half-filled case, shows significant dispersion with 
transfered momentum - the peak disperses outwards in 
energy as we go from the zone center towards the zone 
corner. Finally, we study the RIXS response in the pres- 
ence of finite core-hole broadening. In Fig.[T31[c), we show 
results for a typical broadening (of the order of the Mott 
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FIG. 14: (Color online) Panel (a) : d-hole spectral function, 
j4^(cj), as well as the core-hole spectral function, Aii(lo) with 
(r = 1.0, blue line) and without (F = 0.0, green line) core- 
hole broadening effects in a particle-hole asymmetric large gap 
Mott insulator. For the core-hole spectral function, energy is 
measured with respect to the core-hole energy. Eh, i- e. to ^ 
uj — Eh- Panels (b) and (c), RIXS response as a function of the 
transfered energy, f2, for various fixed transfered momenta, 
X, without [F = 0.0, panel (b)] and with [F = 1.0, panel 
(c)] core-hole broadening effects for incident photon energy 
oji = 3.5. Magnified response in the inset of panel (b) clearly 
shows a two-peak structure and the inset of panel (c) shows 
the quasielastic peak in the blown up far-infrared region. The 
parameters used for this plot are U — 4.0, Qd ~ Qf = 5.0, 
T = 0.1 and n/ = 0.25. 



gap) F = 1.0. The first noticeable feature is the reemer- 
gence of the two-peak structure. This is mainly due to 
huge suppression of the sharp resonating peak in Ah{uj), 
as can be observed in Fig. [TW a). which in effect dras- 
tically reduces the resonating response across the Mott 
gap. The high-energy peak still shows significant disper- 
sion across the entire Brillouin zone and the low-energy 
peak also shows dispersive features. Also, as shown in the 
inset of Fig. fT4lc) , a very weak quasielastic peak similar 
to the half-filled case emerges which is related to the long 
tail in the core-hole spectral function. 



VII. CONCLUSIONS 



In conclusion, we have studied the RIXS response in a 
Mott insulator which is modelled by the Falicov-Kimball 
model. We have considered both the particle-hole sym- 
metric half-filled case as well as the general particle-hole 
asymmetric case. We find that when the incident photon 
energy is lying within the upper Hubbard band, the res- 
onant response shows a two-peak structure arising from 
the intraband (low-energy peak) and interband (high- 
energy peak) relaxation processes (as expected since the 
"gap" is larger but the temperature is the same as be- 
fore). The high-energy peak is found to be much larger 
and sometimes overwhelmingly larger (away from half- 
filling case) than the low-energy peak and shows disper- 
sive features throughout the entire Brillouin zone, while 
the low-energy peak remains more or less nondispersive. 
These distinctive features have already been observed in 
a large class of transition metal K-edge RIXS responses in 
a wide class of oxide materials and have been attributed 
to the nonlocal nature of the Mott gap excitations. 

We also have considered moderately large core-hole 
broadening effects on the RIXS response and we see 
that despite significant change in the RIXS response 
many interesting qualitative features like the two-peak 
structure and the dispersive nature of the high-energy 
peak remains more or less intact. The quasielastic 
feature near the zone center becomes comparable to the 
other two peaks. However, this peak will be completely 
masked by the resolution limited tail of the elastic peak 
(and cannot be observed in any current experiments). 
For the half-filled case, we also have calculated the dy- 
namical structure factor, 5'(q, fi), which is proportional 
to the nonresonant part of the response. 5'(q, ft) is either 
dominated by a very weakly dispersive midinfrared peak 
when the transfered photon momentum, X, is near the 
zone corner or by a narrow peak around = when the 
momentum is close to the zone center and exactly at the 
zone center 5'(q, fl) vanishes for finite frequencies, while 
the narrow peak becomes a delta function peak. 
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